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Abstract 

The effects of trimer continuum resonances are considered in the three-body recom- 
bination rate of a Bose system at finite energies for large and negative two-body 
scattering lengths (a). The thermal average of the rate allows to apply our formula 
to Bose gases at ultra- low temperatures. We found a good quantitative description 
of the experimental three-body recombination length of cesium atoms to deeply 
bound molecules up to 500 nK. Consistent with the experimental data, the increase 
of the temperature moves the resonance peak of the three-body recombination rate 
to lower values of \a\ exhibiting a saturation behavior. 
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Efimov states [1] result from a general phenomenon that occurs in three-boson 
systems, corresponding to the appearance of an increasing number of bound 
states as the two-body binding energy tends to zero. The number is infinite 
when the two-body energy is zero. So, for systems close to an infinite two-body 
scattering length some curious three-boson dynamical effects are expected to 
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appear due to the existence of weakly bound trimers or continuum resonances. 
Presently the universal properties of three-boson systems can be exploited in 
its plenitude in ultracold atomic traps by varying from positive to negative 
values the two-body scattering length (a) near a Feshbach resonance. 

Recent experimental observation of a signature of Efimov physics was reported 
in Ref. [2], considering an ultracold gas of cesium atoms in a trap. Such results 
definitely open a window to study the universal behavior of three-body systems 
with large scattering lengths in a controlled way, and can be considered as 
an experimental verification of correlations among few-body observables with 
short range interactions with respect to the scattering length. 

A giant three-body recombination loss for a large and negative two-body scat- 
tering length, where the three-body state hits the continuum threshold [3], 
was reported in Ref. [2] . For a particular value of the scattering length, char- 
acteristic of their experiment, the three-body state attains zero binding energy 
and produces a huge increase in the measured recombination rate when three 
free-atoms with low kinetic energy collide forming deeply bound dimers. The 
increase of the gas temperature T from 10 nK to 250 nK depletes the recom- 
bination peak moving it to lower values of \a\. This behavior indicates that 
the zero energy bound state evolves to a triatomic continuum resonance [4] , in 
agreement with results obtained for the trajectory of an Efimov state in the 
complex plane, i.e., the zero energy trimer dives into the continuum forming 
a resonance when \a\ is increased [4]. 

Theoretically, the position and width of the triatomic resonance moves to 
larger values as \a\ diminishes, starting from the one that allows a zero energy 
state [4]. This agrees with the observation that as the temperature raises the 
peak moves to smaller \a\. The average energy available for the recombina- 
tion process raises with temperature, consequently the recombination peak 
(roughly at the resonance position) moves toward smaller values of \a\ and the 
intensity diminishes due to the increase of the decay rate to the continuum 
states. The resonance decay to the three-body continuum states competes with 
the formation of the deeply bound dimer. 

The theory for calculating the zero temperature recombination rate [3,5] rep- 
resents fairly well the 10 nK data. The issue now is how to describe quanti- 
tatively the temperature dependence of the recombination rate and the peak 
dislocation toward smaller \a\ by raising T. 

In the present work, we generalize the formula of Braaten and Hammer for the 
recombination rate into deeply bound states. The expression, derived for the 
particular case of zero kinetic energies [3,5] is extended to positive energies. 
We include the effect of continuum trimer resonances for a < in the recombi- 
nation rate. A thermal average of the recombination rate is performed before 
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quantitative comparison with the experimental results of Kraemer et al. [2] for 
the cesium gas at T = 10, 200 and 250 nK. The key point to generalize the rate 
for positive energies, is to write down the energy and width of continuum res- 
onances in terms of a < and a three-body scale. To obtain the detailed form 
of these functions, we solved numerically the renormalized zero-range three- 
body model [4,6,7] in the complex plane to compute the continuum resonance 
energy and width as a function of a and a three-body scale. 

In practice it is needed one three-body information to set the scale of a real 
system. For example, in the experiment of Ref. [2] performed with a 10 nK 
trapped cesium atoms near a Feshbach resonance there is an evidence of the 
virtual formation of a zero energy state, during the three-atom recombination 
process to a deeply bound dimer. The recombination peak was found at the 
particular value of a in their fitting procedure. This information is enough to 
set the three-body scale. 

The three-body physical scale is given by the binding energy (B 3 ) of the 
shallowest trimer just below the zero-energy Efimov state. The trimer binding 
energy B 3 is calculated from the universal relation [4] 



aZ 1 = 0.0297^/-£?3 , (1) 

where h 2 /m = 48.12/m KA 2 , and m is the boson mass number. In the fol- 
lowing general discussion, the three-body scale will be given by the particular 
value of the scattering length — a_, where the three-boson system has an ex- 
cited Efimov state at zero energy. 

The recombination rate, for a finite energy E, can be written as 



ha 4 

L 3 (E) = 3C(a,E) — , (2) 
m 

where C(a, E) at zero energy is given by [5] : 

C(a, 0) = 4590 , ^W^l ■ ^ ' ® 

sin [sq m(Ja|/a_JJ + smn r/_ 



where is a dimensionless quantity that characterizes the resonance width. 
It describes the unknown decay rate of Efimov states into deeply bound dimer 
states plus a free atom. 

Equations (2), with L 3 = L 3 (0), and (3) were used in Ref. [2] to fit the 
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recombination length p 3 at zero energy, as 



2m 



(4) 



The finite energy extension of Eq. (3) is done by the inclusion of the pole 
brought by the continuum resonance. 

With the trimer at zero energy, the S-matrix pole is obtained from 



C(a,a_) = sin[s m(|a|/a_)], (5) 

when a = — a_ (£(a, a_) = 0). Here and in the following, a_ = a_(E = 0) 
to simplify the notation. In a naive extension of the above function to the 
complex energy plane, the corresponding zeros, £(a, a_(£ 1 )) = 0, will give the 
energy and partial width of the continuum resonances related to the scattering 
length a = —a-(E). 

In extending Eq. (3) to non-zero energies, it is assumed that the parameter r}_, 
which characterizes the width of the resonance, varies slowly with the scatter- 
ing energy. This implies that it is not expected an appreciable effect when the 
trimer moves from the threshold to an ultralow energy continuum resonance. 
However, the calculated decay width of the resonances to continuum states 
has a magnitude large enough to be important for the observed results. 

The slow variation of in terms of the energy is expected because it corre- 
sponds to the formation of a deeply bound dimer with a size much smaller 
than \a\. Therefore, the short-range part of the triatomic wave function should 
be responsible for the recombination process, which is deeply inside the poten- 
tial, such that does not strongly depend on the wave function tail, determined 
by the large scattering length. 

The problem now is focused on how to get a_ as a function of energy. Or, how 
to obtain a_ (E) . The resonance energy and partial width scale with a and a_ 
(a_ gives the dependence on the three-body scale). Therefore, we can write 
that 



E = T T{aja_) , (6) 

m at 

where the scaling function T[ajaS) should be universal in the limit of large 
values of a and a_, with respect to the interaction range. Note that the imag- 
inary part of E in Eq.(6) gives the decay width to the continuum states. For 
a = —a- the trimer has zero energy implying that T{— 1) = 0. By expanding 
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the real and imaginary parts of T{z) around z — — 1, we obtain the approxi- 
mate formula 



ma 



h 2 



E 




(7) 



where a and 7 are fitting parameters. This expression is found to reach reason- 
able results when comparing to the numerical results of the model presented 
in Ref. [4]. 

In figure 1 we show the dimensionless products = Re(e) and ej = Im(e) as 
functions of a/a_, for the three-body energy scale derived from the condition 
of an Efimov state at zero energy for a = — a_ . The simple parameterization is 
compared to the calculated values using a renormalized zero-range three-boson 
model [4], with the three-body scale adjusted to allow an Efimov state at zero 
energy. The real part of the resonance energy behaves almost linearly for the 
range of values of a/a_ shown in the figure, while the imaginary part of e 
has an almost quadratic dependence on a/a_. The values of the dimensionless 
constants of Eq. (7), a = 1.01 and 7 = 3.58, are chosen to fit the numerical 
results obtained with the renormalized zero range model of Ref. [4]. 




Fig. 1. Real, 6r (lower frame (b)), and imaginary, ej (upper frame (a)), parts of 
the dimensionless expression for e = E{ma 2 _/h 2 ) as a function of (a/a_). The solid 
lines are obtained from the renormalized zero model [4] and the dashed lines are 
the corresponding analytic results given by Eq. (7) with a = 1.01 and 7 = 3.58. 



One can easily find the function a_(E) by solving Eq.(7) for a given value of 
E, with a = —a-(E): 
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a.(E) 



1 + i 



a 

Y 1 



7 ma_ 
Ai--r—-E 



cr 



h 2 



(8) 



Equation (8) is strictly valid for E being the resonance energy with the imagi- 
nary part related to the partial width to the continuum states, in spite of that 
we perform a naive analytic continuation of Eq. (8) to values of E along the 
real axis. Then, instead of (5), we will have 



((a, a-(E)) = sin 



s In 



(9) 



with 



The function ((a, a -(E)) has a zero at the resonance energy for a = —a -(E) 
when the decay width to the deep states vanishes. This extends the Efimov 
law of appearance of bound states to the case of continuum resonances as a 
function of the scattering length. Scaling a by exp (tt/sq), and considering that 
a and 7 are scale invariants, the zero of ((a-(E), a) gives a resonance energy 
re-scaled by a factor of exp (— 2-n/ sq). 

The recombination rate for finite energy is derived from the original expression 
at E = by introduction of the complex function for ((ci-(E),a), given by 
Eq. (9), assuming that the main effect for the three-body recombination into 
deep dimer states comes from the change of the nearest S-matrix pole due to 
the triatomic continuum resonances, while the rj- is a slowly varying function 
of E. The resulting analytic expression for L 3 (E) in Eq. (2) is given by 



C(a, E) = 4590 ,,, , r , 2 - (11) 

v ' ; \((a-(E),a) +isinhr/_| 2 v ' 



The next step is to perform a thermal average of the recombination rate (L^)t- 
Considering that the bosons are not condensed, we will use a classical thermal 
energy distribution, which allows the separation of the center of mass motion. 
The remaining degrees of freedom in the center of mass kinetic energy E are 
the Jacobi energies of the two independent relative motions, then E = E p + E q , 
where the index p and q are related, respectively, to the relative momentum of 
a pair of atoms, p, and the relative momentum of the third atom to the pair, 
q. The thermal average of the recombination rate calculated classically, with 
P = l/(kT) and k the Boltzmann constant, is given by 
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(L 3 ) T = V J dE v ^E p J dE q yjE~ q L 3 (E P + E q ) e ~^ +E ^ . (12) 



The three-body observables are correlated to the two- and three-body ref- 
erence scales in the Thomas- Efimov limit, where the ratio a/ro (r^ is the 
effective range) goes to infinity (scaling limit). These correlations are gener- 
ally defined in terms of functions of dimensionless quantities named scaling 
functions. Using a_ as a reference, an observable can be written in terms of a 
scaling function [7], such that 



0(E) = (a_)^(e,a/a_) , (13) 

where O is a general observable of the three-body system at energy E with 
the power v giving the proper dimension to the observable. The function 
JF (e, a/a__) converges to a limiting curve when the Thomas-Efimov condition 
is approached [7]. The scaling functions, calculated in the the scaling limit [8], 
are identified with a limit cycle [9]. 

From the above Eq. (13), the thermal average of the scaling function of a given 
observable is 



(0(E)) T = a»_(F(e,a/a-)) T = (a_)"^ 



th 



ma 



h 2 



■kT, a/ a. 



(14) 



which defines a correlation between a and T for a given value of (0)t with 
the reference scale a_. The observed peak position of the recombination rate 
as a function of T is one of such example. 

The parameters of the simple curves for the resonance energy enters through 
Eq. (11) in the calculation of the thermal average of the recombination length, 



2m 



(15) 



In figure 2, the results of the present model are compared to the experimental 
data given in [2]. We allowed a change in the parameter, from 0.06 found 
in [2] to 0.03, as we have also considered a_ to be 0.898a instead of 0.85a - 
Replacing a_ and the cesium mass number m = 133 in Eq. (1) we found 
B 3 = 1.81 mK for the experiment. The experimental data has temperatures 
around 10 nK for one set (full circles), 200 nK (full triangles) and 250 nK 
(open diamonds). Our result for 200 nK is somewhat above the data, while 
for 300 nK approaches the experimental values. The increase of temperature 
makes the resonance peak lower and wider, almost disappearing for 500 nK. 
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Fig. 2. Recombination length {pz)t [as given by Eq. (15)] in a trapped gas of cesium 
atoms, as a function of the scattering length and temperature. The solid curves from 
up to bottom are the theoretical results for T = 10 nK, 100 nK, 200 nK, 300 nK, 
400 nK and 500 nK. The symbols are the experimental results for T = 10 nK (full 
circles), 200 nK (full triangles) and 250 nK (open diamonds) from [2]. 
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Fig. 3. Position of the maximum of the recombination length as a function of the 
temperature. For the experimental data, see comments in [10]. 

The movement of the peak to lower values of \a\ is also seen in the figure. The 
present simplified model of introducing temperature dependence in the three- 
body recombination rate is also compatible with results obtained in Ref. [11]. 

In figure 3, we show the results for the scattering length corresponding to the 
peak position as a function of temperature compared to recent and preliminary 
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experimental results [10]. In the present model, no new parameters are used 
to obtain these results. 

In summary, for a three-boson system, it was derived an extension to positive 
energies of an expression given in [3] for the recombination length, including 
the triatomic continuum resonance pole. The position and width of the tri- 
atomic resonance increase as the modulus of the scattering length is decreased 
near the critical condition for which an Efimov state turns into a continuum 
resonance. 

The thermal average of the observables defines new scaling functions appro- 
priate to ultracold gases. Our approach was tested with experimental data of 
Ref. [2], as shown in Figs. 2 and 3. We found that for the particular thermal 
average of the recombination rate into deep dimer states, the gross features of 
the experimental data for T below 500 nK are reproduced, with parameters 
set at T = 0. The scattering length for the peak position as a function of 
temperature describes nicely the data and the saturation for T above 200 nK 
is well described. 

To conclude, the trajectory of the triatomic continuum resonance is reflected 
in the correlation curve of a and T, which is a direct consequence of dominance 
of just two scales: the two and three body ones which determine the three-body 
observables. Therefore, the thermal scaling law was checked by the experiment 
with cesium atoms which strongly supports the dominance of few scales in the 
physics of ultracold atoms. Such universal thermal scaling law for three-body 
recombination length into deep dimers are expected to be confirmed also in 
other systems. 

In view of the nice results obtained in comparison with cesium experimental 
data, this simplified analytical approach of introducing temperature depen- 
dence in few-body thermal observables for ultracold atomic gases is shown to 
give the basic properties of a more detailed formalism to be developed. 

We thank Dr. R. Grimm and T. Kraemer for useful discussion and for pro- 
viding us preliminary data referred in [10]. This work was partially supported 
by Fundagao de Amparo a Pesquisa do Estado de Sao Paulo and Conselho 
Nacional de Desenvolvimento Cientifico e Tecnologico. 
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